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Abstract 


High levels of theory have been used to compute quartic force fields (QFFs) for the cyclic and 
linear forms of the C .H . + molecular cation, referred to as c-C ,H , + and 1-C 3 H 3 + . Specifically the 
singles and doubles coupled-cluster method that includes a perturbational estimate of connected 
triple excitations, CCSD(T), has been used in conjunction with extrapolation to the one-particle 
basis set limit and corrections for scalar relativity and core correlation have been included. The 
QFFs have been used to compute highly accurate fundamental vibrational frequencies and other 
spectroscopic constants using both vibrational 2 nd -order perturbation theory and variational 
methods to solve the nuclear Schrodinger equation. Agreement between our best computed 
fundamental vibrational frequencies and recent infrared photodissociation experiments is 
reasonable for most bands, but there are a few exceptions. Possible sources for the discrepancies 
are discussed. We determine the energy difference between the cyclic and linear forms of C 3 H 3 + , 
obtaining 27.9 kcal/mol at 0 K, which should be the most reliable available. It is expected that 
the fundamental vibrational frequencies and spectroscopic constants presented here for c-C 3 H 3 + 
and 1-C 3 H 3 + are the most reliable available for the free gas-phase species and it is hoped that 
these will be useful in the assignment of future high-resolution laboratory experiments or 
astronomical observations. 
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1. Introduction 


The two lowest energy forms of C 3 H 3 + are cyclopropenyl cation and propargyl cation, though 
the energy difference between these is large - with the cyclic form being approximately 26 
kcal/mol lower in energy. 1 Cyclopropenyl cation is the smallest aromatic carbocation, which 
explains its stability, but it should be noted that even with this large energy difference, it is 
common to find both isomers when they are produced in the gas-phase, suggesting that they form 
from different mechanisms. There has been considerable experimental and theoretical work on 
both isomers, and we refer the interested reader to Refs. 1-9 and references therein. Here, we 
discuss a few of the earlier works that relate to this study, but first we note that for convenience 
we shall refer to cyclopropenyl cation as c-C 3 H 3 + and propargyl cation as 1-C 3 H 3 + , and when 
referring to both isomers we use C 3 H 3 + . 

Our interest in c-C 3 H 3 + stems from astrochemistry. Cyclopropenylidene, c-C 3 H 2 , which 
possesses a large dipole moment, has been shown to be ubiquitous in the interstellar medium 
(ISM), 10 12 and its main formation pathway has been proposed to be due to dissociative 
recombination of an electron with c-C 3 H 3 + . 13-14 Hence there has been considerable interest in 
detecting c-C 3 H 3 + in the ISM for more than two decades. However, c-C 3 H 3 + possesses D <h 
symmetry, resulting in no permanent dipole moment, and thus is not detectable via microwave 
(rotational) spectroscopy. Furthermore, there is no experimental high-resolution rovibrational 
spectrum available to analyze astronomical observations. In fact, until recently there was no gas- 
phase spectrum of its vibrational frequencies, but instead only matrix isolation spectra or spectra 
from salts. 1516 That changed in 2002 when Dopfer et al 2 4 used infrared photodissociation (IRDP) 
experiments to observe the C-H stretching region of C 3 H 3 + complexed with various ligands. In 
2010, Ricks et al 5 improved upon these experiments by measuring the gas-phase infrared (IR) 
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spectrum of the isomers of C,H 3 + that were associated with one Ar atom. The results of these 
latter two studies are generally consistent with the matrix isolation experiments, though many 
more bands were assigned in the Ricks et al IRPD experiment, going down to approximately 
1100 cm 1 . One inconsistency that was noted by Ricks et al, however, was that the assignment 
for the doubly degenerate CH stretching mode, v 4 (e'), at 3182 cm 1 was about 44 cm 1 higher than 
theory. They attributed this to problems in scaling factors for the theoretical calculations, and 
suggested that further theoretical work was needed. Thus one of the purposes of the present 
study was to provide theoretical predictions of the fundamental vibrational frequencies of both 
C .H , + wherein scaling is not required. 

Determination of an anharmonic force field and fundamental vibrational frequencies has 
actually been reported for c-C 3 H 3 + in two 1989 studies by Lee et al 6 and by Xie and Boggs. 7 Lee 
et al computed a full quartic force field (QFF) at the Hartree-Fock level of theory to determine 
anharmonic corrections via second-order rovibrational perturbation theory which they then 
appended to harmonic frequencies computed at the second-order Mpller-Plesset (MP2) 
pertubation level of theory. They reported an array of spectroscopic constants from their 
anharmonic analysis including anharmonic constants, vibration-rotation interaction constants, 
and quartic and sextic centrifugal distortion constants for c-C 3 H 3 + and its deuterated 
isotopologues. One interesting issue they uncovered was the fact that the standard formula used 
to compute the spectroscopic constants for symmetric top molecules 1719 can be incorrect when 
there is a non-totally symmetric, non-degenerate vibrational mode, as there is for c-C 3 H 3 + . The 
issue was discussed in some detail and Lee et al confirmed the modifications they made to the 
standard symmetric top formula by slightly perturbing the mass of one atom thereby slightly 
breaking the Du, symmetry and forcing the SPECTRO program 20 to use an asymmetric top 
analysis. Xie and Boggs used the MP4(SDQ) level of theory (fourth-order Moller-Plesset 
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perturbation theory including single, double, and quadruple excitations, but excluding the more 
expensive triple excitations) to construct a partial potential energy surface that included diagonal 
force constants through sixth-order, but only second-order off diagonal force constants (in a local 
internal coordinate system). They then used an approximate variational approach to solve the 
nuclear Schrodinger equation. A detailed comparison of the two studies was presented in the 
Lee et al paper. In short, for the levels of theory used, both studies gave reasonable agreement 
with the matrix isolation experiments available at the time and with each other, though it is 
interesting to note that Lee et al obtained 3178 cm" 1 for V 4 (e') while Xie and Boggs obtained 
3149 cm" 1 . 

The higher energy isomer, 1-C3H3 + , has also received considerable attention from both theory 
and experiment. It was observed in both of the IRPD experiments mentioned above, 2 ' 5 and it has 
also been studied recently at high levels of theory by Botschwina and Oswald. 8 Botschwina and 
Oswald used an explicitly correlated method, CCSD(T)-F12x (x=a,b), 21-22 which is based on the 
singles and doubles coupled-cluster method that includes a perturbational correction for triple 
excitations, denoted CCSD(T). 23 They computed a five dimensional potential energy function, 
involving the totally symmetric modes, which included up through sixth-order diagonal constants 
and up through fourth-order off diagonal force constants, and solved the nuclear Schrodinger 
equation variationally. The five dimensions included the symmetric CH 2 stretch, the free 
acetylenic CH stretch, the CH 2 scissor mode, and the two CC stretching modes (a correction was 
applied to the antisymmetric CH 2 fundamental to account for the neglect of anharmonic 
coupling). Botschwina and Oswald found generally good agreement with the IRPD experiment 
of Ricks et al. 

More recently Botschwina and Oswald used the CCSD(T)-F12x (x=a,b) levels of theory to 
examine the equilibrium structure and harmonic vibrational frequencies of Ar complexes of both 
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c-C 3 H 3 + and 1-C 3 H 3 + . 9 They located three distinct minima for Ar complexes of both isomers, 
though one of the C s minima for 1-C 3 H 3 + was noticeably lower in energy than the other two 
minima. For the c-C 3 H 3 + isomer, all three minima were relatively close in energy. The results of 
this latter study reaffirmed Botschwina and Oswald’s assignment of the 3238 cm 1 band in Ricks 
et al’s IRPD experiments to the acetylenic CH stretch vibration in 1-C 3 H 3 + . 

The accuracy of state-of-the art ab initio predictions for fundamental vibrational frequencies 
has improved considerably in the last twenty years, and it is common now to determine 
fundamental vibrational frequencies to within a few wavenumbers (cm 1 ) of high-resolution 
experiments (for example, see Refs. 24-27). To this end, theoretical spectroscopists predicted 
more than a decade ago 28 that state-of-the art ab initio predictions were becoming reliable 
enough that it should be possible to assign an astronomical spectrum using only ab initio 
predictions and without high-resolution laboratory experimental data for difficult species, such as 
small molecular anions and cations. This has now occurred only a few years ago when 
Cernicharo et al 29 reported detecting the small molecular anion C 5 N in the C-rich star 
IRC+10216 and based their assignment on the ab initio calculations of Aoki 30 and Botschwina 
and Oswald. 31 

Thus our goal in the present study is to compute highly accurate QFFs for the c-C 3 H 3 + and 1- 
C 3 H 3 + molecules, and to predict their rovibrational spectroscopic constants to very high accuracy. 
These data are of interest to astronomers now more than ever given that the Herschel Space 
Observatory is in operation and collecting high-resolution data, the NASA Stratospheric 
Observatory for Infrared Astronomy (SOFIA) has begun its series of initial science flights, and 
the James Webb Space Telescope (JWST), often referred to as the replacement to the Hubble 
Space Telescope, will launch later in this decade. Further, the Atacama Large Millimeter Array 
(ALMA) is set to start early science operations in late 2011. Some of the instruments for these 
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telescopes operate at longer wavelengths and thus are not useful for rovibrational spectroscopy 
of c-C ,H , + , but they may be able to detect deuterium and l3 C isotopologues due to their increased 
sensitivity, which is especially interesting for isotopologues of c-C 3 H 3 + . The parent isotopologue 
has no permanent dipole moment, hence it cannot be observed via rotational spectroscopy, but 
isotopologues that do not retain D }h symmetry will exhibit a small permanent dipole moment 
since the center of nuclear charge and center of mass of the molecule will no longer be the same 
(and the molecule rotates about its center of mass). For 1-C 3 H 3 + , the lowest energy vibration may 
be within range for instruments on all of the above telescopes, but since 1-C 3 H 3 + possesses a 
permanent dipole moment, the parent isotopologue as well as all deuterium and 13 C 
isotopologues may be detectable. However, a full spectroscopic analysis of all deuterium and 
l3 C isotopologues of c-C 3 H 3 + and 1-C 3 H 3 + is beyond the scope of the present study, and will be 
reported separately. 32 

The Theoretical Approach is described in the next section, followed by Results and 
Discussion. Our Conclusions are presented in the final section. 

2. Theoretical Approach 

A. Details of the Electronic Structure Methods 

We first describe details of the electronic structure calculations, including details of the 
corrections that have been included. In general, we follow the approach we have developed in 
recent years 33,34 in which we extrapolate CCSD(T) energies to the one -particle basis set limit, 35 
followed by addition of corrections for scalar relativity 36 and core correlation. The valence 
CCSD(T) calculations were performed in conjunction with Dunning’s correlation consistent 
basis sets. 37 We will denote the cc-pVXZ (X=T, Q, or 5) basis sets as TZ, QZ, or 5Z. A 
correction for scalar relativity is evaluated at the CCSD(T)/TZ level of theory using the Douglas- 
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Kroll approximation. 36 As we pointed out previously, 33 the scalar relativity integrals lose 
precision when going beyond the TZ basis set, which is problematic for computing QFFs. Core 
correlation was included as a correction by performing CCSD(T) calculations, with and without 
the core correlated, using the Martin-Taylor basis set designed for this purpose. 38 

QFFs have been determined according to the prescription described previously. 33,34 For both c- 
C 3 H 3 + and 1-C 3 H 3 + , a reference geometry was determined at the CCSD(T)/5Z level of theory with 
corrections for core correlation and scalar relativity taken into account. A grid of displacement 
geometries centered on this reference structure (and based on the symmetry internal coordinates 
discussed later) was then used for all calculations. The number of unique geometries was 1961 
and 2479 for c-C 3 H 3 + and 1-C 3 H 3 + , respectively. As indicated, CCSD(T) energies are 
extrapolated to the one-particle basis set limit using a three-point formula that experience has 
shown to be reliable. 33,34,39 The scalar relativity and core-correlation corrections are added, and 
the energies are used in a least squares fit of a QFF for each molecule. For c-C 3 H 3 + , the 1961 
unique energies were augmented to a redundant set of 3837 energies that was used to fit 460 
unique coefficients in the QFF. The sum of the squared residuals was 1.31 x 10 17 a.u. 2 . For 1- 
C 3 H 3 + , the 2479 unique energies were augmented to a redundant set of 4565 energies that was 
used to fit 572 unique coefficients in the QFF (sum of the squared residuals = 2.61 x 10 17 a.u. 2 ). 
The final QFF for each molecule was then obtained by an analytical transformation to the exact 
minimum (i.e., to where the gradient terms are exactly zero). We note that it has been shown that 
some molecules with C-C multiple bonds exhibit erratic behavior for bending frequencies, 
including molecules like acetylene, ethylene, and benzene, 39 43 although cyclopropenylidene does 
not show this behavior. 44 The problem is associated with ensuring that the one-particle basis set 
is properly balanced with respect to saturation in the lower angular momentum functions (i.e., s 
and p functions) and inclusion of higher angular momentum functions. Though we report only 
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our best QFF here for both c-C,H , + and 1-C,H , + , we have examined in detail several QFFs for 
both isomers and found that neither suffers from this issue. All electronic structure calculations 
were performed with the MOLPRO 2006.1 program. 45 

B. Details of the Coordinate Systems and the Vibrational Methods 

For both molecules, the QFFs were determined in symmetry internal coordinates. For 1-C 3 H 3 + , 

we use the following definition of symmetry internal coordinates: 

•Si(Aj) = R 4 (C - H); 

S 2 (A 1 ) = (R la + RJ/^2- 

5 3 (Aj) = (R 2 [C = C] + R 3 [C = C])/V 2; 

5 4 (Aj) = {R 2 [C = C] - R 3 [C = CJ)/v/2; 

S s (A l ) = O l (AH a CH h )- 
S 6 (B l ) = (R la -R lb )/^; 

5 7 (Bj) = UN\(C -C-C- rout)\ 

5 8 (B,) = LIN\(C -C-H- r ou t); 

S 9 ( BQ = (LHJCC - AH h CC)/^2; 

5 10 (B 2 ) = LIN\(C -C-C- rin); 

S n (B 2 ) = LIN\{C -C-H- r m ); 

S l2 (B 2 ) = r(H a -C-C-H b ) 

where the simple internal coordinates for 1-C'3H 3 + are given in Fig. 1 . Note that Sj and Sh are in- 
plane bending modes and Sw and Sn are out-of plane bending modes. The definition of LINl(a- 
b-c-d) and %{a-b-c-d) are taken from the INTDER program: 46 

LIN\(a - b - c - d) = sin -1 [e,/ • (ebc x eta )] ; 

T (a-b-c-d) = sm~\e ba -{e c b x e C d) / (sin (p ahc - sin (p bcd )] 

where the e are unit vectors defined as e a b = et, — e a . The reference vectors r,„ and r 0 „t were 
defined as follows: 
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y in — R 2 X ( R 1 a X R 1 b ) 
y out = Rj X y in = R 2 X (R 2 X ( R\a X R\b )) 

Fig. 1 also gives the planar equilibrium and ground state (GS) vibrationally averaged structures 
for both c-C 3 H 3 + and 1-C 3 H 3 + , as well as the GS vibrationally averaged rotational constants. The 
vibrationally averaged quantities are the “position average,” i.e., r z , computed from 2 nd -order 
perturbation theory. Symmetry relationships for the quadratic, cubic, and quartic force constants 
are given later. 

For c-C 3 H 3 + , we use the exact same definition of symmetry internal coordinates as given by 
Lee et al. 6 They are repeated here for convenience: 

5) (flj ) = (i?! + R 2 + R 3 ) / V3 , 

S 2 (a l ) = (q +r, + r 3 )/V 3, 

S 3 (o 2 ) = (cq + a 3 + a 5 - a 2 - a 4 - a 6 ) / V6 , 

S 4 c( e ) = Oi ~r 2 -r 3 )/V 6, 

5 5a (e) = ( ■2p i -0 2 -P 3 )/S, 

S 6a (e ) = (aj + a 2 + a 4 + a 5 - 2 a 3 - 2 a 6 ) / VT2 , 

5 7 (a 2 ") = (y, +y 2 +/ 3 )/V3, 

5 8a (e") = (2 7l -y 2 -y 3 )/V6, 

= (^2 -r 3 )/j2, 

S 5h (e) = (P 2 -P 3 )/S, 

S bh (e) = (a 2 +a 5 -a \ -a 4 )/ 2, 

S sb (e'') = ( Y2 -y 3 )/y[2, 

where the simple internal coordinates R, r, a, /i are the bond lengths and bond angles defined in 
Fig. 2, and y refers to the out of plane bending angle for a given C-H bond with respect to the 
plane defined by the three C atoms (see Fig. 2). 

Fundamental vibrational frequencies were computed using either a vibrational variational 
method (VAR) or second-order perturbation theory (PT). 17 19 The MULTIMODE program 47 was 
used for the VAR calculations, while the SPECTRO program 20 was used for the PT calculations, 
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and for computing other spectroscopic constants. For 1-C.H, + , the QFF was analytically 
transformed into a simple-internal, Morse-cosine coordinate system for the vibrational 
variational calculations. The benefits of using Morse coordinates for the stretches when using a 
QFF in vibrational variational calculations can be traced back to Meyer et al 48 and Carter and 
Handy 49 more than 20 years ago. Later in 1994, Dateo et al 50 first defined the Morse (3 parameter 
solely on the computed force constants (i.e. (3=-F iii /(3F ii )) instead of optimizing it with respect to 
experimental data. We follow this (3 definition, which requires that the transformed diagonal 
cubic force constant for the stretch vanish. For c-C 3 H 3 + , it is necessary to use a symmetry 
adapted Morse-cosine coordinate system (for the stretches and bends) because it is a ringed 
system. In addition, instead of the out-of-plane coordinates (S 7 , S 8a , S 8b ), we use the sine of these 
coordinates. Rather than transforming the symmetry internal coordinate QFF into the symmetry 
adapted Morse-cosine-sine coordinate system, we refit the QFF. We note that transformation of 
the QFF into a Morse-cosine coordinate system is important (and in the case of c-C 3 H 3 + , a 
Morse-cosine-sine coordinate system), otherwise some fundamental vibrational frequencies, in 
particular stretching frequencies, can be too high by tens of cm 1 . As discussed in Ref. 50, the 
Morse-cosine coordinate system serves to build in the correct limiting behavior for the potential 
function. 

3. Results and Discussion 

A, Equilibrium Structures and Harmonic Frequencies 

The equilibrium structure, equilibrium rotational constants, and harmonic frequencies for c- 
C 3 H 3 + and 1-C 3 H 3 + are presented in Tables 1 and 2, respectively. The equilibrium structures 
obtained in this work have somewhat shorter bond distances relative to those published 
recently 89 using CCSD(T*)-F12a, and the HCH angle determined here for 1-C 3 H 3 + is about 1° 
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larger. We note that the longer C-C bond in 1-C 3 H 3 + , R 2 (see Fig. 1), is intermediate between a 
typical double and single bond length, but it is much closer to that in ethylene rather than ethane. 
Conversely, the corresponding harmonic frequency o) 4 is intermediate between a typical double 
bond and a single bond, but its value is closer to that for ethane rather than ethylene. Previous 
authors have referred to this C-C bond as a “single bond,” so we adopt that terminology here, but 
we note that it is intermediate between a single and double bond, which also means that there 
will not be free rotation of the terminal CH 2 group about this bond. The shorter C-C bond is very 
much like a typical triple bond both in its bond length and harmonic frequency. 

The harmonic frequencies obtained here for 1-C 3 H 3 + are in reasonable agreement with those 
given by Botschwina and Oswald 8 , though our stretching harmonic frequencies are generally a 
few cm 1 higher, consistent with the shorter bond lengths obtained in the present work. We note 
that we use a different convention for the symmetry labeling of the modes relative to Refs. 5 and 
8 - essentially the B[ and B 2 labels are reversed. That is, following the convention that Herzberg 
used for C 2v planar molecules, B[ is used for in-plane antisymmetric modes and B 2 is reserved for 
out-of-plane vibrations, and that is the convention adopted here. For c-C 3 H 3 + , the harmonic 
frequencies given in Table 1 are in reasonable agreement with the MP2 values from Ref. 6 and 
the MP4(SDQ) values of Ref. 7, given the differences in levels of theory. 

As indicated previously, the energy difference between c-C 3 H 3 + and 1-C 3 H 3 + is about 26 
kcal/mol, but the levels of theory used in the present study should yield a much more definitive 
value. The energy that we obtain at the minimum is -115.7647467662 E h and -115.717377491 
E h for c-C 3 H 3 + and 1-C 3 H 3 + , respectively. The electronic energy difference is thus 10,396.4 cm 1 . 
The anharmonic zero-point energies given by 2 nd -order PT are 9841.5 cm 1 and 9208.0 cm' 1 , 
which includes the E 0 term (the polyatomic equivalent of the a 0 Dunham coefficient for 
diatomics). 51 The corresponding MULTIMODE zero-point energies are 9823.7 cm 1 and 9189.1 
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cm 1 , for c-C.H, + and 1-C,H, + , respectively. The 0 K energy difference we obtain is 27.9 
kcal/mol, which is slightly larger than what has been reported experimentally. 1 While the 
experimental value is not at 0 K, it also has a few kcal/mol uncertainty - see Ref. 1 and 
references therein for more details. The 27.9 kcal/mol OK value obtained in the present work 
should be the most reliable available. 

B. Fundamental Vibrational Frequencies and Spectroscopic Constants 

The GS vibrationally averaged structure and rotational constants, and the fundamental 
vibrational frequencies obtained for c-C ,H , + in the present work are presented in Table 3. Other 
spectroscopic constants obtained from 2 nd order perturbation theory are presented in Table 4 
(anharmonic constants) and Table 5 (vibration-rotation interaction constants, and quartic and 
sextic centrifugal distortion constants). For the variational calculations, we used four mode 
coupling and five mode coupling in order to demonstrate the convergence. Comparison of the 
fundamental vibrational frequencies for the two columns labeled VCI 4MR and VCI 5MR shows 
that indeed there is excellent convergence, with the largest difference being only 1.5 cm 1 for v 7 . 
4MR/5MR refer to the number of modes coupled in the potential energy expansion formula, 
while 4-mode coupling was adopted in all Coriolis integrations. For most vibrational modes, the 
difference is 1.1 cm 1 or less, and for three of the modes the difference is less than 1 cm 1 . Based 
on these comparisons and experience, we would estimate that the variational fundamentals are 
converged to better than 1.0 cm 1 for the VCI 5MR values. Comparison of the VCI 5MR results 
with the fundamentals obtained from 2 nd -order perturbation theory shows good agreement with 
the two approaches, with the differences being consistent to what we usually find for tightly 
bound molecules that do not possess a large amplitude motion. Specifically, the largest 
difference is 3.8 cm 1 for v 2 , but this mode is affected by a Fermi type 1 resonance with 2v 7 . For 
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most of the vibrational modes the difference is less than 3 cm' 1 , again showing that 2 nd -order 
perturbation theory is a good approximation for solving the nuclear Schrodinger equation for a 
tightly bound molecule. 

Our best results should be the VCI 5MR fundamentals. Comparison of these to the available 
experimental data shows reasonable agreement for the matrix isolation values where we might 
expect differences of up to 20 cm 1 or so due to a matrix shift. In fact, the largest differences 
between the matrix isolation results and our VCI 5MR values are only 9.6 and 12.0 cm 1 for v 3 
and v 8 , which are both determined indirectly (see Refs. 15 and 16 for details). Comparison of the 
VCI 5MR results to the IRPD values from Ricks et al 5 shows very good agreement for v 5 , but for 
v 4 we obtain a value that is 47.2 cm 1 lower than their assignment at 3182 cm 1 . Thus our best 
estimate for v 4 is consistent with previous theory and calls into question their assignment. 
Further, the value we obtain for v 4 , 3134.8 cm 1 , is more consistent with the assignment from 
Dopfer et al. 2 4 One of the reasons we performed the variational calculations in the present study 
was to be certain that we had v 4 described properly since our 2 nd -order perturbation theory results 
did not agree with the assignment from Ricks et al. Given the levels of theory used in the present 
study and the fact that we have ruled out any possible resonance issues in solving the vibration 
problem, we can definitively conclude that the 3182 cm 1 band observed by Ricks et al is either 
not representative of the free gas-phase spectrum for c-C 3 H 3 + or it is due to a different vibrational 
mode or species. We have examined the variational Cl results for possible combination bands or 
overtones in the variational calculations that might explain the band at 3182 cm 1 , but none 
appear for either c-C 3 H 3 + or 1-C 3 H 3 + , at least not within 10 cm' 1 . There is a doubly degenerate 
band involving three quanta, 2v 6 +v 5 , that is very close to 3182 cm' 1 , but this seems unlikely. It 
may be that the band observed at 3182 cm 1 is shifted somewhat due to complexation with the Ar 
atom, or it may be due to a different species. 
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Given the levels of theory used in the present study, the spectroscopic constants presented for 
c-C 3 H 3 + in Tables 4 and 5 should be highly accurate, and it is hoped these will be useful in the 
future assignment of high-resolution rovibrational spectra from either laboratory experiments or 
astronomical observations. 

The vibrationally averaged structure and rotational constants, and the fundamental vibrational 
frequencies obtained for 1-C 3 H 3 + in the present work are presented in Table 6. These are 
“position averaged” values (i.e., r z ) computed with 2 nd -order perturbation theory. Other 
spectroscopic constants obtained from 2 nd order perturbation theory are presented in Table 7 
(anharmonic constants) and Table 8 (vibration-rotation interaction constants, and quartic and 
sextic centrifugal distortion constants). Comparison of the VCI 4MR and 5MR results contained 
in Table 6 shows that the variational calculations are converged to better than 1 cm 1 , similar to 
the situation for c-C 3 H 3 + . In fact, the largest difference is only 0.8 cm 1 for v 7 . Comparison of 
the VCI 5MR and 2 nd -order perturbation theory fundamental vibrational frequencies for 1-C 3 H 3 + 
shows reasonable agreement, though not as good as found for c-C 3 H 3 + . The largest differences 
occur for the C-H stretches v t and v 6 , being 10.3 and 9.1 cm' 1 , respectively. We note that v 6 is 
involved with a significant Fermi type 2 resonance with v 3 + v 7 , though the difference between 
the two components of the resonance is fairly consistent between 2 nd -order perturbation theory 
and VCI 5MR (35.7 versus 34.0 cm 1 ). Differences between 2 nd -order perturbation theory and 
VCI 5MR for the other fundamental vibrational frequencies are more in line with the differences 
we found for c-C 3 H 3 + . Interestingly, the agreement between 2 nd -order perturbation theory and 
VCI 5MR for the 2v 4 overtone and the v 12 + v 4 combination band is not nearly as good, which is 
expected as one moves into the realm of less pure states and stronger coupling. 

Agreement between our best VCI 5MR fundamental vibrational frequencies and those 
obtained by Botschwina and Oswald 8 using the CCSD(T*)F-12a/VTZ-F12 level of theory is 
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modest. The largest differences occur for Vj (15 cm 1 ) and v 4 (27 cm 1 ). There are many 
possible sources for these differences, with the most significant probably being the use of an 
approximate (T) contribution by Botschwina and Oswald, their neglect of core correlation, and 
their neglect of most coupling to non-totally symmetric vibrational degrees of freedom in solving 
the variational nuclear Schrodinger equation. 

Comparison of our best VCI 5MR results with the experiments of Ricks et al shows reasonable 
agreement for most of the assignments with a few exceptions. Our best value for v, agrees very 
well, confirming the conclusion by Botschwina and Oswald 9 that the 3238 cm 1 band observed by 
Ricks et al is essentially a free acetylenic C-H stretch. The agreement for v 2 , v 3 , and v 4 is also 
very good, with differences all less than about 10 cm 1 . Agreement for the overtone band 2v 4 is 
reasonable, being about twice the difference for the v 4 fundamental, and agreement for v 6 is also 
reasonable as the assignment by Ricks et al falls between the two components of the Fermi type 
2 resonance between v 6 and v 3 + v 7 . Agreement for v 5 , v 10 , and the combination band v I2 + v 4 is 
more modest, however. Botschwina and Oswald 8 have already questioned the reliability of the 
assignment for the totally symmetric mode v 5 . Given that the lowest energy structure found for 
1-C 3 H 3 + • Ar has the Ar atom out of plane and over the C-C single bond, and that this structure is 
quite a bit lower in energy than the other minima, 9 it seems plausible that the C-C single bond 
stretch v 5 would be significantly impacted in the complex. This same reasoning could be applied 
to v l0 , which is described as a CFF out-of plane wag, and to the combination band v 12 + v 4 since 
v 12 is an out-of plane bending mode that involves the CCC backbone. Thus, the discrepancies 
found between the assignments of Ricks et al and our VCI 5MR results can reasonably be 
attributed to shifts in the vibrational frequencies as a result of complexation for 1-C 3 H 3 + . 

o 

For /-C 3 H 3 , the effects of scalar relativistic corrections are small: ~ -5E-5 A on the C-H bond 
lengths and — 3E-4 A on the CC bonds; less than or equal to 0.4 cm on harmonic frequencies; 
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and -6 to 0 cm 1 for the vibrational fundamentals. As expected from previous experience, core- 

o o 

correlation effects are much larger on the geometry (-1.3E-3 A for the CH bonds and -3E-3 A for 
the CC bonds) and harmonic frequencies (+2 to +7 cm 4 ), but only slightly larger for the 
fundamentals (+2 to +8 cm' 1 , except +10.7 cm 1 for v 10 ). For c-C 3 H 3 + , the effects of scalar 

o o 

relativity corrections are smaller than found for Z-C 3 H 3 : -2.1E-4 A for R cc ; -5.5E-5 A for r CH ; 
0.0-0. 5 cm 1 for harmonic frequencies; and -1 to +2 cm 1 for fundamentals (except +8 cm 1 for 

o o 

v 3 ). However, core-correlation effects are larger: -5.4E-3 A for R cc ; -1.4E-3 A for r CH ; +2 to +7 
cm 1 for harmonic frequencies; and +3 to +9 cm 1 for fundamentals (except +0.3 cm 1 for v 6 ). 
More details are available upon request. 

Given the levels of theory used in the present study, the spectroscopic constants presented for 
1-C 3 H 3 + in Tables 7 and 8 should be highly accurate, and it is hoped these will be useful in the 
future assignment of high-resolution rovibrational spectra from either laboratory experiments or 
astronomical observations. 


C. Quartic Force Fields 

For completeness, the best QFFs computed in this work are given in Tables 9 through 12. 
Specifically, Table 9 contains the quadratic and cubic force constants and Table 10 the quartic 
force constants for c-C 3 H 3 + . Table 1 1 contains the quadratic and cubic force constants and Table 
12 the quartic force constants for 1-C 3 H 3 + . These are given in symmetry internal coordinates and 
symmetry relationships between the force constants are given in the tables. The force constants 
presented are based on the following quartic expansion: 


^ i.j i.j.k 


1 


H 

24 
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where the summations are unrestricted. We note that for the force constant labels for 1-C .H . + , we 
have labeled modes 10, 11, and 12 as x, y, and z, respectively, in order to avoid confusion. So, 
for example, the diagonal quadratic force constant for mode 10 is given as F xx . 

4. Conclusions 

Accurate CCSD(T) QFFs have been computed for the c-C 3 H 3 + and 1-C 3 H 3 + molecular cations. 
Extrapolation to the one-particle basis set limit has been included as well as corrections for scalar 
relativity and core correlation. Anharmonic spectroscopic constants have been determined from 
2 nd -order perturbation theory and fundamental vibrational frequencies have been evaluated from 
2 nd -order perturbation theory and from variational calculations. Agreement between 2 nd -order 
perturbation theory and variational Cl calculations for the fundamental vibrational frequencies is 
very good. Agreement between our computed fundamental vibrational frequencies and recent 
IRPD experiments is good with a few exceptions. The 3182 cm 1 band assigned in one recent 
IRPD experiment 5 to v 4 for c-C 3 H 3 + does not agree well with our calculations where we obtain 
3134.8 cm 1 (VCI 5MR), but our value does agree well with other recent IRPD experiments. 2 4 
We have examined our VCI calculations for possible combination and overtone bands, but can 
find only one band that could be a reasonable match and this requires three quanta. Hence we 
conclude that the 3182 cm 1 band may be a combination or overtone band that is perturbed 
somewhat by the presence of the Ar atom, or it may be due to a different species. For 1-C 3 H 3 + , 
agreement between the IRPD experiments and our variational calculations for the fundamental 
vibrational frequencies (plus one overtone and one combination band) is reasonable except for 
v 5 , v 10 , and the combination band v 12 + v 4 . However, a recent ab initio study 9 that explored the 
potential energy surface of 1-C 3 H 3 + interacting with an Ar atom shows one minimum quite a bit 
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lower than the others, and this structure would likely exhibit perturbations to v 5 , v 10 , and the 
combination band v l2 + v 4 , so a plausible explanation for these discrepancies is given. 

The fundamental vibrational frequencies and spectroscopic constants presented here for c- 
C 3 H 3 + and 1-C 3 H 3 + should be the most reliable available for the free gas-phase species and it is 
hoped that they will be useful in the assignment of future high-resolution laboratory experiments 
or astronomical observations. Finally, we compute what should be the most reliable energy 
difference between the c-C 3 H 3 + and 1-C 3 H 3 + , obtaining a value of 27.9 kcal/mol at 0 K. 
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Figures and Tables 


Figure 1. Equilibrium structures, and vibrationally averaged structures (italic numbers) and 
rotational constants of c-C 3 H 3 + and 1-C3H3 + determined from our best QFF. See text for details. 
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Fig. 2. Internal coordinate definitions used for the c-C3H3 + cation. 
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Table 1. Equilibrium structure (A / deg), rotational constants (cm 1 ), and harmonic frequencies 
(cm' 1 ) for c-C 3 H 3 + , determined from our best QFF (CCSD(T)/cc-pV(T,Q,5)Z extrapolation + 
core-correlation + scalar relativistic corrections). See text for more details. 



Geometry 


Harmonic Vibrational Frequencies 



Parameters 

Mode 

Symmetry and Description 

Freq 

R( C-C) 

1.3590363 

0)1 

a{ (breathing, C-H str - C-C str) 

3312.3 

KC-H) 

1.0777461 

0) 2 

a{ (breathing, C-H str + C-C str) 

1650.8 

A e = B e 

1.03260 

0)3 

a{ (in-plane internal torsion) 

1058.2 

C e 

0.51630 

0) 4 

e (in-plane deformation) 

3265.1 

Z.C-C-H 

150.0 

0) 5 

e' (in-plane wagging) 

1326.7 

zc - c-c 

60.0 

0) 6 

e' (in-plane scissoring) 

946.0 



0) 7 

a{ (symmetric out-of-plane bends) 

764.5 



0) 8 

e" (asymmetric out-of-plane bends) 

1023.4 


Table 2 . Equilibrium structure (A / deg), rotational constants (cm 1 ), and harmonic frequencies 
(cm' 1 ) for /-C 3 H 3 + , determined from our best QFF (CCSD(T)/cc-pV(T,Q, 5 )Z extrapolation + 
core-correlation + scalar relativistic corrections). See text and Fig. 1 for more details and the 
definition of R1-R4 and 0 \ . 


Geometry 

Parameters 

Mode 

Harmonic Vibrational Frequencies 
Symmetry and Description 

Freq 

Ri 

1.085506 

0)1 

a\ (C-H str) 

3367.15 

Ri 

1.345571 

0)2 

a\ (CH 2 sym str) 

3122.73 

Ri 

1.227265 

0)3 

a\ (C=C str) 

2131.15 

R4 

1.072298 

0) 4 

a\ (CH 2 bending) 

1483.44 


120.355 

0) 5 

ci\ (C-C str) 

1137.60 

Z.C-C-H 

180.0 

0) 6 

b\ (CH 2 anti-sym str) 

3232.42 

zc-c-c 

180.0 

0) 7 

b\ CH 2 group rocking 

1039.76 



0)g 

b\ (CCH in-plane linear bend) 

631.00 

A e 

9.53209 

0) 9 

b\ (CCC in-plane linear bend) 

292.51 

B e 

0.32329 

O)l0 

b 2 (CH 2 out-of-plane wag) 

1120.50 

C e 

0.31269 

0)11 

b 2 (CCH out-of-plane linear bend) 

882.00 



0)12 

b 2 (CCC out-of-plane linear bend) 

254.71 


Table 3. Vibrationally averaged structure (A / deg), rotational constants (cm 1 ), and vibrational 
fundamental frequencies (cm' 1 ) for c-C 3 H 3 + , determined from our best QFF. 2 nd -order 
Perturbation Theory (2 nd PT) and vibrational variational Cl results, where nMR represents the 
highest mode coupling level in the potential term expansions, are presented and compared to 
experiment. See text for more details. 


Zero-Point Structure and Rotational Constants 


Rz(C-C) 

1.3658234 


X 

o 

II 

to 

o 

1.02610 



Rz(C-H) 

1.0806204 


Co 

0.51178 




Anharmonic Vibrational Analysis 



Mode 

2 nd -PT 

VCI 

4MR 

VCI 

5MR 

Exp c 

Exp d 

Exp 

Vl(Ol') 

3176.6 
->3171. T 

3174.8 

3175.4 

3183 



v 2 (ai') 

1618.3 
->1620. 7 b 

1622.0 

1622.1 

1626 



v 3 (a 2 ') 

1040.3 

1039.3 

1040.6 

(1031) 


3125±4 e 

Me') 

3131.7 

3133.8 

3134.8 

3138 

3182 

3 1 29 f 







3130 s 

Me') 

1299.6 

1295.9 

1296.2 

1290 

1293 


Me') 

924.2 

925.9 

927.0 

927 



v 7 (a 2 ") 

756.6 

755.6 

757.1 

758 



Me") 

1004.5 

1000.9 

1002.0 

(990) 



ZPE 

9841.5 

9833.2 

9834.0 





a Fermi type 1 resonance with 2v 2 
b Fermi type 1 resonance with 2v 7 


c Refs. 15 and 16. 
d Ref. 5. 
e Ref. 2. 
f Ref. 3. 
s Ref. 4. 


Table 4. The x st and g„anharmonic constant matrices for c-C 3 H 3 + determined from our best 
QFF. All values are in cm' 1 . See text for more details. 


Mode 

1 

2 

3 

4 

5 

6 

7 

8 

1 

-18.820 








2 

-2.749 

-3.596 







3 

-6.862 

-3.610 

-1.030 






4 

-75.767 

-2.134 

-5.084 

-28.685 





5 

1.634 

-12.513 

-6.429 

0.943 

-6.667 




6 

-7.881 

-2.876 

-1.420 

-6.993 

-3.418 

-1.647 



7 

-7.703 

-3.452 

0.727 

-6.373 

-0.076 

-1.681 

-1.811 


8 

-7.001 

-2.323 

3.143 

-5.993 

1.398 

1.159 

-1.541 

0.077 


Mode 

4 

5 

6 

8 

4 

9.532 




5 

-0.684 

3.384 



6 

0.064 

-0.361 

0.762 


8 

-0.073 

-0.287 

-0.182 

0.578 


Table 5. Vibration-rotation interaction constants and quartic and sextic centrifugal distorsion 
constants for c-C 3 H 3 + . See text for more details. 


Vib-rot Constant (MHz) 
Mode a B a c 

Distortion Constant 
(MHz) 

(Hz) 


Watson S reduction 
(MHz) 

(Hz) 

1 

92.1 

45.9 

X aaaa 

-0.291 

OCaa 

0.183 

Dj 

0.073 

Hj 

0.259 

2 

90.0 

44.9 

"C bbbb 

-0.291 

‘hbbb 

0.335 

D jk 

-0.122 

H jk 

-1.119 

3 

-13. 8 a 

16.6 

X cccc 

-0.023 

*hccc 

0.001 

D k 

0.055 

Hkj 

1.466 

4 

85.7 

42.2 

X aabb 

-0.291 

‘I’aab 

0.956 

di 

0.000 

H k 

-0.605 

5 

-15.2 

42.1 

X aacc 

-0.047 

‘habb 

-0.179 

d2 

0.000 

h, 

0.000 

6 

137. 8 a 

21.7 

X bbcc 

-0.047 

*haac 

-0.171 



h 2 

0.000 

7 

-109. 8 a 

-10.9 



*hacc 

0.002 



h 3 

-0.038 

8 

23. 7 a 

-26.4 



^bcc 

0.002 










‘hbbc 

-0.171 










‘habc 

-0.342 






a Modes for which Coriolis resonance has been taken into account. 


Table 6. Vibrationally averaged structure (A / deg), rotational constants (cm 1 ), and vibrational 
fundamental frequencies (cm 4 ) for 1 -C 3 H 3 + , determined from our best QFF. 2 nd -order 
Perturbation Theory (2 nd PT) and vibrational variational Cl results, where nMR represents the 
highest mode coupling level in the potential term expansions, are presented and compared to 
experiment. See text for more details and Fig. 1 for coordinate definition. 


Zero-Point Structure and Rotational Constants 

Ri 

1.095770 


6 x 

120.284 

R2 

1.353993 


Aq 

9.40357 

S3 

1.236271 


So 

0.31956 

R4 

1.068242 


Co 

0.30861 


Anhannonic Vibrational Analysis 



2 nd PT 

Variational Cl 

Exp d 



4MR 

5MR 


Vi (Ai) 

3228.7 

3238.9 

3239.0 

3238 
3 1 39 e 

V 2 a (dO 

2997.0 

2999.2 

2998.7 

3004 

V3 04 1) 

2084.0 

2082.2 

2082.2 

2077 

v 4 (A 1 ) 

1429.8 

1433.7 

1434.4 

1445 

v 5 b (d,) 

1128.5 

1131.9 

1131.8 

1222 

v 6 c (B 1) 

3061.9 

3070.8 

3071.0 

3093 

3097.6 

3104.5 

3105.0 

3 1 1 3 f 

V7 (Si) 

998.0 

999.8 

1000.6 


v 8 (Si) 

598.0 

607.8 

607.7 


v 9 (SO 

294.0 

294.2 

294.8 


V10 (S 2 ) 

1054.6 

1057.9 

1058.1 

1111 

V11 (S 2 ) 

859.7 

861.8 

861.9 


Vi 2 (S 2 ) 

249.3 

251.8 

251.7 


2v 4 

2836.9 

2856.8 

2857.8 

2878 

Vi 2 +V 4 

1676.4 

1695.2 

1695.5 

1755 

V5+V3 

3193.4 

3202.0 

3201.8 

3 1 9 l/3243 f 
3184/3238® 


a Fermi resonance Type 1 with 2v 4 
b Fermi resonance Type 2 with v u +vi 2 
c Fermi resonance Type 2 with V 3 +V 7 

d Ref. 5. 
e Ref. 2. 
f Ref. 3. 
s Ref. 4. 


Table 7 . The matrix of anharmonic constants x sl for /-C3H3 + , determined from our best QFF. All 
values are in cm’ 1 . See text for more details. 


Mode 

1 

2 

3 

4 

5 

6 

7 

8 

9 10 

11 

1 

- 54.889 










2 

- 0.146 

- 27.831 









3 

- 5.779 

- 1.751 

- 8.478 








4 

- 1.230 

- 31.736 

- 3.781 

- 2.102 







5 

1.635 

- 2.248 

- 15.680 

- 3.110 

- 0.881 






6 

- 0.335 

- 113.563 

- 0.901 

- 21.757 

- 1.314 

- 32.215 





7 

- 1.719 

- 11.150 

- 2.502 

- 17.450 

- 2.595 

- 11.681 

- 3.844 




8 

- 20.642 

- 2.628 

- 11.108 

- 6.735 

- 2.675 

- 2.796 

- 8.808 

- 5.238 



9 

- 2.301 

- 0.575 

- 8.408 

- 2.241 

3.420 

- 0.498 

- 1.749 

- 2.636 

- 0.679 


10 

- 3.339 

- 10.399 

- 6.230 

- 5.004 

- 3.113 

- 20.36 

- 4.794 

- 16.356 

- 2.625 - 9.210 


11 

- 20.948 

- 1.829 

0.992 

- 2.515 

- 3.018 

- 1.529 

- 3.121 

27.704 

1.906 - 13.415 

- 6.040 

12 

- 2.493 

- 1.797 

- 5.193 

- 3.296 

7.240 

- 1.374 

- 2.602 

1.710 

21.446 - 9.393 

- 4.744 


Table 8. Vibration-rotation interaction constants and quartic and sextic centrifugal distorsion 
constants for 1-C3H3 + . See text for more details. 


Mode 

Vib-rot constants (MHz) 
a A a B a c 

Distortion Constant 
(MHz) (Hz) 

Watson S reduction 
(MHz) (Hz) 

1 

-21.277 

25.909 

24.454 

^ aaaa 

-85.377 

^aaa 

4809.680 

Dj 

0.003 

Hj 

-0.001 

2 

4879.106 

6.846 

11.688 

bbbb 

-0.012 

d^bbb 

0.000 

D jk 

0.479 

H jk 

4.988 

3 

151.243 

55.273 

51.836 

t’ 

L cccc 

-0.010 

^CCC 

0.000 

D k 

20.862 

Hkj 

-2267.423 

4 

-3188.803 

3.755 a 

11.156 

^ aabb 

-1.931 

^aab 

-742.911 

d ' i 

- 0.000 

H k 

7072.116 

5 

225.856 

23.754 a 

24.150 a 

'E aacc 

-0.007 

^abb 

5.534 

(%2 

- 0.000 

hi 

- 0.000 

6 

3052.329 

8.691 

9.786 

^ bbcc 

-0.011 

^aac 

-1515.360 



h 2 

0.001 

7 

-3310. 921 a 

-8.912 

6.57 l a 



^acc 

-0.059 



h 3 

0.000 

8 

336.274 a 

-12.094 

-0.485 



^bcc 

0.000 





9 

-2468. 800 a 

-54.337 

-24.327 a 



‘t’bbc 

-0.001 





10 

3490. 356 a 

0.795 a 

-5.110 



^abc 

5.962 





11 

403.182 a 

2.727 

-6.288 









12 

1278. 499 a 

-42.247 a 

-62.886 










a Modes for which Coriolis resonance has been taken into account. 


Table 9. Complete set of non-zero quadratic and cubic force constants for c-C 3 H 3 + in a symmetry 
internal coordinate system. See text for more details. Units of force constants are mdyn/A n .rad ra 
appropriate for an energy unit of mdyn A ( 1 mdyn A = 1 aJ). 


Quadratic and Cubic Force Constants 


F„ 

7.217432 

F44i=F4b4bl 

-0.1448 

F655 = -F 65b5b = -F 6b5b5 

-0.4172 

F 2 i 

-0.098532 

F 442 = F 4b4b2 

-19.1998 

F65b3 = -F6b53 

0.3817 

F22 

5.798770 

F 444~-F 44b4b 

-13.5720 

F 65b4b = F 6b54b 

-0.1160 

F33 

0.262156 

F 541 = F 5b4bl 

0.2160 

F661 = F6b6bl 

-0.1679 

F 44 = F 4b4b 

5.800512 

F 542 = F 5b4b2 

-0.1399 

F 662 = F 6b6b2 

-0.0729 

F 54 = F 5b4b 

-0.009810 

F 544 = -F 54b4b = -F 5b44b 

-0.0774 

F 664 = -F 6b64b = -F 6b6b4 

0.0495 

F55 = F5b5b 

5.316006 

F 54b3 = "F 5b43 

-0.1386 

F665 = -F 6b65b = -F 6b6b5 

-0.2352 

F 64 = F 6b4b 

0.017754 

F551 = F5b5bl 

-13.225 

F 666 F 6b6b6 

0.0959 

F 65 = F 6b5b 

0.298226 

F 552 = F 5b5b2 

-0.2622 

F771 

-0.2385 

F 66 F 6b6b 

0.414829 

F 554 = -F5b54b = -F 5b5b4 

-0.1110 

F772 

-0.0827 

F77 

0.371862 

F555 =_ F5b5b5 

-18.9961 

F 874 = F 8b74b 

-0.0760 

F88 = Fgb8b 

0.490412 

F641 = F6b4bl 

-0.0146 

F 875 = F 8b75b 

0.3281 

Fill 

-22.5615 

F 642 = F 6b4b2 

-0.0009 

F 876 = F 8b76b 

0.2836 

F21I 

-0.1640 

F 644 = -F 64b4b = -F 6b4b4 

0.0073 

F 881 = F 8b8bl 

-0.4132 

F22I 

0.1680 

F 64b3 = "F 6b43 

-0.0602 

F 882 = F 8b8b2 

-0.1193 

F222 

-19.2136 

F651 = F6b5bl 

-0.1032 

F 884 =- F 8b84b = -F 8b8b4 

-0.0776 

F33I 

-0.2589 

F 652 = F 6b5b2 

-0.1790 

F 885 =- F 8b85b = -F 8b8b5 

0.3106 

F332 

-0.0596 

F 654 = F 6b5b4 

0.1160 

F 886 =- F 8b86b = -F 8b8b6 

0.2418 


PH PP 


Table 10. Complete set of non-zero quartic force constants for c-C 3 H 3 + in a symmetry internal 
coordinate system. See text for more details. Units of force constants are mdyn/A n .rad ra 
appropriate for an energy unit of mdyn A ( 1 mdyn A = 1 aJ). 


Quartic Force Constants 


Fill! 

59.44 

F 7744“ F 7 7 4b4b 

F2111 

0.20 

F 7754“ F 775b4b 

F2211 

-0.45 

F 7755 = F 775b5b 

F2221 

-0.35 

F 7764 = F 776b4b 

F2222 

56.66 

F 7765 = F 776b5b 

F3311 

0.02 

F 7766“F 776b6b 

F3321 

0.11 

F 8741“F 8b74bl 

F3322 

-0.00 

F 8742 = F 8b74b2 

F3333 

0.04 

F8751 = Fgb75bl 

F7711 

-0.11 

F 8752 = F 8b75b2 

F7721 

0.15 

F 8761“F 8b76bl 

F7722 

-0.10 

F 8762“F 8b76b2 

F7733 

0.67 

F 881 1 = F 8b8bl 1 

F 774b3 

-0.02 

F8821 = Fsb8b21 

F775b3 

0.01 

F 8822 = F 8b8b22 

F776b3 

0.04 

F 8833“F 8b8b33 

F7777 

0.09 

F 8844 = F 8b8b4b4b 

F44n=F4b4bll 

-0.37 

F 884b4b = F 8b8b44 

F 442 1 = F4b4b2 1 

-0.42 

F 8854“F 8b8b5b4b 

F 4422“F 4b4b22 

56.74 

F 8855 = F 8b8b5b5b 

F 4433~F 4b4b33 

-0.00 

F 885b4b = F 8b8b54 

F5411 = F 5 b4bll 

-0.26 

F 885b5b = F 8b8b55 

F 542 1 = F 5b4b2 1 

0.18 

F 8864“F 8b8b6b4b 

F 5422“F 5b4b22 

0.35 

F 8865“F 8b8b6b5b 

F 5433 = F 5b4b33 

-0.02 

F 8866 = F 8b8b6b6b 

F551 i=F 5 b5bl 1 

20.59 

F 886b4b = F 8b8b64 

F5521 = F5b5b21 

0.50 

F 886b5b = F 8b8b65 

F 5522 = F 5b5b22 

-0.34 

F 886b6b = F 8b8b66 

F 5533“F 5b5b33 

-0.14 

F 8877 = F 8b8b77 

F 5544™ F 5b5b44 

-0.33 

F444i=-F 4b4b41 

F 554b4b = F 5b5b4b4b 

-0.34 

F 4442 =- F 4b4b42 

F6411— F6b4bll 

0.05 

F 4b443“"F 4b4b4b3 

F 642 1 “F6b4b2 1 

0.02 

F 54b31“"F 5b43 1 

F 6422 = F 6b4b22 

0.05 

F 54b32 =- F 5b432 

F 6433“F 6b4b33 

0.01 

F5551 = ‘F 5b5b5 1 

F6511 = F6b5bll 

-0.22 

F5552—F 5b5b52 

F6521 = F6b5b21 

0.03 

F 5b553“"F 5b5b5b3 

F 6522“F 6b5b22 

-0.16 

F64b31 = -F 6b431 

F 6533 = F 6b5b33 

0.04 

F 64b32“"F 6b432 

F 6544 = F 6b5b4b4b 

-0.10 

F65b31 = -F 6b531 

F 654b4b = F 6b5b44 

-0.13 

F65b32 = -F 6b532 

F 6554 = F 6b5b5b4b 

-0.13 

F6661 = -F 6b6b61 

F 65b5b4 = F 6b554b 

0.05 

F 6662 =- F 6b6b62 

F6611 = F6b6bll 

0.03 

F 6b663“"F 6b6b6b3 


-0.11 

F 5554 = F 5b5b5b4b = 3 F 5b5b54 = 3 F 5b554b 

-0.50 

-0.04 

F 6444 = F 6b4b4b4b = 3 F 64b4b4 = 3 F 6b4b44 

0.05 

-0.02 

F 6555 = F 6b5b5b5b = 3 F 65b5b5 = 3 F 6b5b55 

0.55 

0.03 

F 6664 = F 6b6b6b4b = 3 F 6b6b64 = 3 F 6b664b 

0.19 

0.03 

F 6665 = F 6b6b6b5b = 3 F 6b6b65 = 3 F 6b665b 

0.19 

0.94 

F544 i=-F 54b4b 1 = "F 5b4b4 1 

0.09 

0.15 

F 5442 = -F 54b4b2 = 'F 5b4b42 

0.21 

-0.09 

F 54b43 = F 5b443 = -F 5b4b4b3 

0.04 

-0.18 

F554 1 = -F 5b5b4 1 = -F 5b54b 1 

0.14 

-0.10 

F 5542 = -F 5b5b42 = 'F 5b542 

0.02 

-0.11 

F 554b3 = F 5b543 = ‘F 5b5b4b3 

0.28 

-0.08 

F644 1 = -F 64b4b 1 = -F 6b4b4 1 

-0.01 

-0.08 

F 6442 = -F 64b4b2 = -F 6b4b42 

0.03 

0.26 

F 64b43~F 6b443 = ‘F 6b4b4b3 

0.05 

-0.11 

F 655 1~"F 65b5bl = "F 6b5b5 1 

-0.10 

0.62 

F 6552 = -F 65b5b2 = 'F 6b5b52 

0.26 

-0.14 

F 65b53 = F 6b553 =_ F 6b5b5b3 

0.24 

-0.01 

F664 1 = -F 6b6b4 1 = "F 6b64b 1 

-0.04 

0.03 

F 6642 = -F 6b6b42 = "F 6b64b2 

0.11 

-0.45 

F 664b3 = F 6b643 = ‘F 6b6b4b3 

-0.02 

-0.11 

F665 1 = 'F 6b6b5 1 = -F 6b65b 1 

0.01 

0.55 

F 6652 = -F 6b6b52 = "F 6b65b2 

0.06 

-0.03 

F 665b3 = F 6b653 =_ F 6b6b5b3 

0.08 

-0.20 

F 8744““F 874b4b = "F 8b74b4 

-0.05 

0.43 

F 8755 =- F 875b5b = ”F 8b75b5 

-0.76 

-0.00 

F 8766 =- F 876b6b = ”F 8b76b6 

-0.68 

0.22 

F 884l”“F 8b84bl = "F 8b8b41 

0.12 

1.39 

F 8842 =- F 8b84b2 = "F 8b8b42 

-0.06 

0.47 

F 884b3 = F 8b843 = ‘F 8b8b4b3 

-0.03 

-0.30 

F 8851~"F 8b85bl = "F 8b8b5 1 

-0.41 

40.12 

F 8852~"F 8b85b2 = "F 8b8b52 

-0.40 

0.00 

F 885b3 = F 8b853 = ‘F 8b8b5b3 

-0.22 

0.14 

F 8861~“F 8b86bl = "F 8b8b61 

-0.02 

-0.08 

F 8862~"F 8b86b2 = "F 8b8b62 

-0.09 

38.85 

F 886b3~F 8b863 = ‘F 8b8b6b3 

-0.53 

0.06 

F654 1 = -F65b4b 1 = -F 6b54b 1 = ‘F 6b5b4 1 

-0.06 

-0.18 

F 6542 = -F 65b4b2 = "F 6b54b2 = ‘F 6b5b42 

0.07 

0.10 

F 654b3 = F 65b43 = F 6b543 = ‘F 6b5b4b3 

0.08 

-0.03 

F 8754 =- F 875b4b = ”F 8b754b = ‘F 8b75b4 

0.02 

-0.06 

F 8764 =- F 876b4b = ”F 8b764b = ‘F 8b76b4 

0.03 

-0.12 

F 8765 =- F 876b5b = -F 8765b = ‘F 8b76b5 

-0.68 

0.06 

12 Non-Symmetry-Unique Constants: 

-0.08 

F7654— (F6644-F6655)/2 

0.00 

0.01 

F 9854 = (F 8844"F 8855)/2 

0.01 


F6621 = F6b6b21 

0.24 

F 874b3“"F 8b743 

-0.04 

F9876 = (F 8 866-F 8 877)/2 

0.00 

F 6622~F 6b6b22 

-0.11 

F 875b3“"F 8b753 

-0.11 

F 65b4b4 = F 6b54b4 = (F 6544‘F 654b4b)/2 

0.02 

F 6633“F 6b6b33 

0.09 

F 876b3“"F 8b763 

-0.79 

F 65b54b = F 6b5b54 = (F 6554"F 65b5b4)/2 

-0.09 

F 6644 = F 6b6b4b4b 

-0.09 

F 8887 =- F 8b8b87 

0.57 

F 6b65b4 = F 6b654b = (F 6654‘F 665b4b)/2 

0.06 

F 664b4b = F 6b6b44 

-0.12 

F 4444 = F 4b4b4b4b = 3 F 4b4b44 

85.18 

F 8b84b4 = (F 8 844-F 8 84b4b)/2 

-0.07 

F 6654“ F 6b6b5b4b 

-0.02 

F 5555 = F 5b5b5b5b = 3 F 5b5b55 

92.62 

F 8b85b5 = (F 8855“F 885b5b)/2 

-0.50 

F 6655 = F 6b6b5b5b 

0.07 

F 6666 F 6b6b6b6b 3F 6b6b66 

-0.16 

F 8b86b6 = (F 8 866-F 8 86b6b)/2 

-0.48 

F 665b4b = F 6b6b54 

-0.14 

F 8888~F 8b8b8b8b = 3F 8 b8b88 

0.43 

F 8b854b = F 8b85b4 = (F8854-F 8 85b4b)/2 

0.07 

F 665b5b = F 6b6b55 

0.06 

F 5444~F 5b4b4b4b = 3 F 5b4b44 = 3 F 54b4b4 

0.55 

F 8b864b = F 8b86b4 = (F8864-F 8 86b4b)/2 
F 8b865b = F 8b86b5 = (F 8 865-F 8 86b5b)/2 

-0.01 

-0.21 


Table 11. Complete set of non-zero quadratic and cubic force constants for 1-C3H3 + in a 
symmetry internal coordinate system. See text for more details. Units of force constants are 
mdyn/A n .rad m appropriate for an energy unit of mdyn A (1 mdyn A = 1 aJ). x/y/z represent the 
10 th , 1 1 th , and 12 th internal coordinates, respectively. 


F i 

6.029788 

F 

A XX 

0.167448 

F444 

15.2311 

F 765 

-0.0249 

F965 

0.0962 

Fyx5 

-0.0131 

F 2 i 

0.005712 

F 

± yx 

0.000441 

F511 

0.0002 

F771 

-0.0363 

F971 

0.0092 

Fyyl 

-0.1937 

F 2 2 

5.540738 

F 

1 yy 

0.385229 

F521 

0.0047 

F772 

0.0272 

F972 

0.0408 

Fyy2 

-0.0120 

F 3 i 

-0.119775 

F 

A zx 

0.000222 

F 5 22 

-0.1527 

F773 

-0.9361 

F973 

0.1397 

Fyy 3 

-0.4253 

f 32 

0.072319 

F 

1 zy 

0.015043 

F531 

-0.0050 

F774 

0.1170 

F974 

-0.0580 

Fyy 4 

0.1286 

F33 

11.885630 

Fzz 

0.238988 

F 5 32 

0.1349 

F 775 

0.0325 

F975 

-0.0546 

Fyy5 

0.0064 

F 4 i 

0.144706 

Fin 

-34.5269 

F533 

0.0086 

F 8 61 

0.0002 

F981 

0.0010 

Fzxl 

-0.0026 

F 42 

0.074156 

F211 

-0.0010 

F 5 41 

-0.0043 

F 8 62 

-0.0015 

F982 

0.0095 

Fzx2 

-0.0063 

F43 

-2.970851 

F 22 i 

-0.0055 

F 5 42 

0.1028 

F 8 63 

0.0073 

F983 

-0.0528 

Fzx3 

-0.0961 

F44 

10.906522 

f 222 

-22.3363 

F543 

0.1047 

F 8 64 

0.0128 

F984 

0.0503 

F zx4 

0.0867 

F51 

0.001714 

F311 

0.1236 

F544 

0.2391 

F 8 65 

-0.0053 

F 985 

-0.0134 

Fzx5 

0.0082 

F 52 

0.065407 

F 32 i 

0.0056 

F55I 

0.0024 

F 87 l 

-0.0002 

F99I 

0.0086 

Fzyl 

-0.0032 

F53 

-0.231773 

f 322 

0.0869 

f 552 

-0.3928 

F 872 

-0.0025 

F992 

-0.2028 

Fzy2 

-0.0007 

F54 

-0.173960 

F331 

-0.0602 

F553 

-0.1270 

F 873 

0.1851 

F993 

-0.1762 

Fzy3 

0.0332 

F55 

0.682966 

F 332 

-0.1223 

F554 

-0.0076 

F 874 

-0.2203 

F994 

0.0320 

Fzy4 

-0.0053 

F66 

5.533706 

F333 

-48.6285 

F555 

-0.2658 

F 875 

0.0052 

F995 

0.3241 

Fzy5 

0.0068 

F76 

-0.002072 

F 4 11 

-0.1359 

F66i 

-0.0068 

F 881 

-0.1404 

Fxxl 

-0.0281 

Fzz, 

-0.0031 

F77 

0.283223 

F421 

-0.0008 

F 66 2 

-22.4743 

F 882 

-0.0192 

F X x2 

-0.0185 

F zz 2 

-0.0768 

F 8 6 

-0.004390 

f 422 

0.1801 

F 663 

0.1739 

F 883 

-0.6080 

F X x3 

-0.3280 

Fzz3 

-0.2289 

F 8 7 

0.109569 

F431 

0.0147 

F 664 

0.3239 

F 884 

0.7960 

F xx4 

-0.0273 

Fzz4 

-0.0271 

F 8 8 

0.227801 

F432 

-0.2336 

F 665 

0.2737 

F 885 

-0.0047 

F X x5 

0.0040 

F zz 5 

0.4596 

F96 

0.109515 

F433 

13.1266 

F76I 

-0.0003 

F96I 

-0.0011 

Fyxl 

0.0022 



F97 

-0.067911 

F441 

0.0753 

F762 

-0.0023 

F962 

0.0428 

Fyx2 

0.0098 



F9 8 

-0.023718 

F442 

-0.5212 

F763 

0.0911 

F963 

-0.2622 

Fy x 3 

0.1577 



F99 

0.514243 

F443 

-45.3965 

F764 

0.0117 

F964 

-0.3053 

Fy x 4 

0.0810 




<L> <D 

O ^ 

a 

*H 

,o 


X 

0) 


a 

<D 

c/3 
<D 

*H 

Oh 
<1> 
J-H 

<D -N 

*5 

S o 

(U '-J 
'4—* C3 
cn 

, " H 

2 111 

.5 c 

"d 

d t3 

8 a 

o ^ 

13 ^ 
£< 
a p 
•£ •£ 
>, s 

+2 c+H 

<1> o 

E +3 

a a 


03 60 

. d 

C <u 
•'"' d 
f u 

§ 

J-H 

,o 


SC 

co 

U 


d <u 

.2 

C/3 C 

P ^ 

E ^ 

-2 fe 

S 3 & 
O 03 


o 

g 


CD 


13 

o 

03 

> 

$-4 

$-< 


£ 

C3 

«< 

*4-» 

o 

a/ 

O 

• 

4 -j » 

$H 

£ 

Gh 

C/3 

<D 

03 

T3 

J-H 

2 

C 

rv 

a" 

£ 

C/3 

<D 

o 

<D 

H 

-4—* 

03 

Vh 

<d 

c3 

.£ 

N 

C/3 

-4 — i 

*3 

g 

o 

£ 

§ 
-4— » 
C/3 

S-i 

O 

O 

o 

<H4 

£ 


o 

O 

13 


o 

g 

a> 

<D 

<D 

C/3 

O 

4-* 

a> 

-4-» 

CD 

J-h 

£ 

.£ 

-c 

"Hh 

c 

<+H 

o 

CN 

£ 

o 

C/3 
-4— » 


U 

'S 

£ 

o3 


(N 


a> 

p a> o 
H t3 


00 

O 

O 

t> 

r- 

>■> 

N 

P-H 

in 

o 

NO 

OO 

>N 

N 

PH 

oro- i8 N 

p 

o 

oo 

oo 

>> 

N 

tin 

CN 

O 

O 

NO 

ON 

>> 

N 

[Jh 

IT3 

O 

O 

l> 

On 

N 

Ph 

O 

oo 

ON 

N 

Ph 

CN 

O 

C) 

ON 

O'. 

>> 

N 

Ph 

O 

O 

X 

X 

N 

Ph 

CN 

O 

O 

X 

Pn 

N 

Ph 

OO'O 

CN 

O 

N 

N 

Ph 

CO 

O 

o 

CN 

N 

N 

Ph 

CN 

O 

CN 

n 

N 

N 

Ph 

(N 

p 

O 

CO 

N 

N 

Ph 

oro Z£ZZ d 

CO 

CO 

p 

CO 

CO 

N 

N 

Ph 

O 

O 

N 

N 

Ph 

p 

o 

CN 

N 

N 

Ph 

^sO 

p 

p 

CO 

N 

N 

Ph 

cn 

p 

p 

N" 

N" 

N 

N 

Ph 

N" 

O 

O 

m 

N 

N 

Ph 

p 

CN 

in 

N 

N 

Ph 

CO 

p 

co 

•n 

N 

N 

Ph 

OS 

p 

p 

N- 

m 

N 

N 

Ph 

CN 

OC 

o 

in 

in 

N 

N 

Ph 

in 

CN 

c> 

NO 

NO 

N 

N 

Ph 

O 

O 

NO 

N 

N 

Ph 

p 

p 

N 

N 

Ph 

CN 

^1" 

o 

(N 

r- 

m 

m 


OO 

i— H 

(N 


^1" 

in 

© 

© 

H 

in 

in 

© 

© 

N" 

in 

<N 

r- 

CN 

© 

© 

i 

o 

o 

o 

© 

© 

© 

© 

i— ' 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

i— ' 

CN 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 


CN 


CN 

m 


CN 

m 

N" 


CN 

CO 


in 

NO 

NO 


NO 

i> 

oo 

NO 


oo 

ON 

X 



CN 











in 

in 

in 

m 

in 



l> 

OO 

oc 

00 





X 





X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 





N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

N 

Pd 

PH 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Ph 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

— 

co 

co 

OO 

© 


© 

OO 

^H 1 

<N 

© 

© 

r- 

^1" 

© 


r- 

© 

in 


© 

i— H 

CO 

© 

<N 

© 

N" 

m 

r- 

o 

i—i 

o 

1—1 

© 

© 

m 

© 

© 

CO 

CN 

© 

i— ' 

CO 

© 

CN 

CO 

© 

co 

1— ' 

© 

1 

r- 

© 

© 

© 

© 

i— < 

i— ' 

o 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 


in 

NO 

NO 

l> 

NO 

o 

oo 

NO 


oo 

ON 

X 



CN 


CN 

co 


CN 

co 



CN 

CO 


in 

NO 

IT) 

m 



l> 

00 

oo 

oo 





X 











in 

in 

in 

in 

in 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 







>> 

>> 




>> 






>n 

>N 

>N 

>N 




>3 

X. 







>3 

>N 


>3 









>N 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

co 

i— H 

^1" 

r- 

H 

<N 

(N 

© 

© 

H 

Of- 

t" 

<N 

^1" 

in 

© 

OO 

<N 

co 

H 

r- 

OC 


H 

© 

© 


co 

r- 

© 

o 

co 

(N 

© 

m 

© 

© 

I/O 

© 

© 

i— ' 

© 

© 

© 

© 

© 

CO 

(N 

© 


© 

N" 

© 

© 

CN 

CN 

CN 

in 

o 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 


© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

l> 

OO 

ON 



CN 


CN 

co 


CN 

CO 



CN 

co 


in 

NO 

NO 

l> 

NO 


OO 

NO 


oo 

ON 

X 

ON 

ON 

ON 


CN 

CN 

m 

CO 

CO 





in 

in 

in 

in 

m 

NO 

l> 

l> 

oo 

oo 

oo 

ON 

ON 

ON 

ON 

X 

ON 

ON 

ON 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

ON 

ON 

On 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

pH 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

pH 

Pd 

Pd 

Pd 

pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

CN 

i— H 

co 

H 

© 

m 

(N 

m 

i— H 

<N 



^1" 

© 

H 

H 

N" 

i— H 

H 

<N 

^1" 

© 


N" 


© 

co 

r- 

in 

o 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

CN 

© 

CO 

© 

CN 

CO 

© 

CN 

in 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 


CN 


CN 

m 


CN 

co 



CN 

co 


in 

NO 

NO 


NO 

t"* 

oo 



CN 


CN 

CO 


CN 

CO 

CN 

CN 

co 

m 

m 




m 

in 

in 

in 

in 

NO 

t"* 

r- 

OO 

OO 

00 


CN 

CN 

CO 

CO 

CO 




oo 

OO 

oo 

oo 

oo 

oo 

OO 

oo 

oo 

oo 

oo 

oo 

oo 

oo 

OO 

OO 

oo 

OO 

oo 

oo 

ON 

On 

On 

ON 

ON 

On 

On 

ON 

On 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

Ph 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Ph 

Ph 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

i— H 


<N 

(N 

© 

© 


<N 

© 

H 

CO 


© 

© 

CO 

© 

<N 

<N 

in 

i" 

H 

CN 

CO 

l-~ 

H 

© 

(N 

© 

o 

CN 

(N 

© 

i— ' 

© 

i— ' 

i— ' 

© 

© 

CO 

© 

© 

© 

© 

© 

CN 

in 

© 

© 

© 

1 

i—i 

CO 

© 

© 

© 

© 

i— ' 

i— ' 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

in 

NO 

NO 

o 

NO 

l> 

OO 



CN 


CN 

co 


CN 

co 



CN 

co 


in 

NO 



CN 


CN 

m 

in 

NO 

o 


OO 

OO 

oo 


CN 

CN 

CO 

CO 

CO 





in 

in 

m 

m 

in 

NO 


CN 

CN 

CO 

CO 

co 

oo 

OO 

oo 

oo 

OO 

oo 

oo 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 



r- 

t"* 

l> 

i> 

oo 

OO 

oo 

oo 

OO 

oo 

oo 

On 

ON 

ON 

On 

ON 

ON 

ON 

On 

On 

ON 

ON 

On 

ON 

ON 

On 

ON 

ON 

ON 

On 

On 

ON 

ON 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

<N 

NO 

(N 

i— H 

i— H 

rn 

m 


l/o 

r-- 

CN 

CO 


© 

© 

r- 

CO 


(N 

^1" 


N" 


CO 

<N 

in 

© 

CN 

O 

o 

O 

© 

© 

© 

© 

© 

© 

© 

© 

i— ' 

© 

© 

© 

© 

© 

© 

© 

© 

CN 

© 

© 

© 

© 

© 

© 

© 

CN 

1— ' 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 


CN 

co 



CN 

m 


in 

NO 



CN 


CN 

co 


CN 

co 



CN 

CO 


m 

NO 

NO 

i> 






in 

m 

in 

in 

m 

NO 


CN 

CN 

co 

CO 

CO 


N" 


N" 

in 

in 

in 

in 

in 

NO 

l> 

o 


NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

t"* 

r- 

l> 

i> 

r- 

t> 

l> 

l> 

c- 

l> 


o 

t"* 



t"* 

r- 

i> 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

oo 

oo 

OO 

oo 

oo 

oo 

OO 

OO 

oo 

OO 

oo 

oo 

oo 

oo 

oo 

OO 

oo 

oo 

oo 

PH 

Pd 

Pd 

Pd 

pH 

Pd 

Pd 

Pd 

Pd 

Ph 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Ph 

Pd 

Pd 

Ph 

pd 

Pd 

pd 

Pd 

Pd 

i"- 

^1" 

<N 

<N 

© 

© 

© 

CO 


^H 1 

CO 

© 

(N 

CO 

1— 4 


1 

m 

1 

y —{ 1 

«n 

OC 

N" 

CO 

CO 

© 

© 

N" 

© 

ON 

t-- 

© 

© 


© 

© 

CO 

© 

© 

i—i 

1— ' 

© 

© 

1 

© 

© 

© 

CN 

© 

CN 

© 

© 

© 

© 

1— ' 

r- 

OC 

© 

© 

00 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

in 

NO 



CN 


CN 

CO 


CN 

co 



CN 

co 


in 

NO 



CN 


CN 

co 


CN 

co 

N" 


in 

NO 


CN 

CN 

m 

m 

CO 





in 

in 

in 

in 

in 

NO 


CN 

CN 

CO 

CO 

CO 





in 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

r- 

l> 


i> 

l> 


i> 


i> 

l> 


NO 

NO 

C- 

l> 

C-" 

o 

t"* 


o 

t"* 

l> 

c- 

I> 

r- 

l> 

l> 

l> 

c- 

t"* 

l> 

t"* 

i> 

f" 


i> 

l> 

i> 

t"* 


Pd 

Pd 

Pd 

Pd 

Pd 

pp 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

pp 

Pd 

Pd 

Pd 

Pd 

Ph 

Pd 

pp 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

H 

CN 

IT) 

© 

H 

OO 

<N 

IT) 

© 

(N 

OO 

*T) 

© 

r- 

(N 

© 

OO 

m 

© 

<N 

CN 

(N 

© 

m 

© 

© 

© 

oc 

CN 

O 

o 

m 

© 

© 

© 

© 

© 

1— ' 

© 

© 

© 

r- 

i— ' 

© 

© 

© 

i— « 

CN 

© 

i— « 

m 

OC 

© 

CN 

i— ' 

CN 

© 

CN 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 


CN 

m 



CN 


CN 

CO 


CN 

CO 




CN 


CN 

CO 


CN 

CO 



CN 

CO 


in 


co 

co 

m 


CN 

CN 

m 

CO 

CO 


N" 

N" 

■st 


CN 

CN 

CO 

CO 

CO 

N - 

4f 

N" 

N" 

in 

in 

in 

in 

m 


co 

co 

m 











in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

in 

NO 

in 

m 

m 

m 

in 

in 

m 

in 

in 

in 

in 

in 

in 

m 

in 

in 

m 

m 

m 

in 

in 

m 

m 

m 

m 

in 

in 

m 

NO 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

O 



i— H 

IT3 

© 

m 

CO 

I/O 

© 

© 

^1" 

© 

<N 

N" 


i— H 

© 

r- 

© 

cn 

CO 

© 

N" 

© 

CN 

CN 

in 

© 

NO 

o 

1— ' 

© 


m 

© 

© 

© 

IT) 

© 


CO 

CN 

00 

CO 

© 

© 

© 

cn 

© 

N" 

CO 

© 

in 

in 

© 

© 

© 


© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

CO 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

© 

r- 




r- 

CN 




CN 



CN 


CN 

CO 




CN 



CN 


CN 

p 

co 



CN 





CN 

CN 



CN 

CN 


CN 

CN 

co 

CO 

CO 



CN 

CN 


CN 

CN 

CO 

CO 

CO 


CN 

CN 

CO 



CN 

CN 

CN 


CN 

CN 

CN 

co 

CO 

co 

CO 

CO 

CO 


CN 

CN 

CN 

co 

co 

CO 

CO 

CO 

CO 

N- 



^ t 


CN 

CN 

CN 

CN 

m 

m 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

N" 


"si- 

N" 


N" 

N" 

N" 

N" 


N" 




Pd 

Pd 

Ph 

Pd 

Ph 

Pd 

Pd 

Pd 

Pd 

Ph 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Ph 

Pd 

Pd 

Pd 

Ph 

Pd 

Ph 

Pd 

Pd 

Pd 

Pd 

Ph 


oo<nooo<n^o<nooon 

ooooooooooooo 


^ot^oo^ot^oot^ x x >> *: 

oo oo oo c-* O'* x n 

NNNNNNNNNNNN 

NNNNNNNNNNNN 

|X,tL, pHfefetUPUfefefeP-lP-l 


£U 


nr - r - inn ^’-^ rH ^ r-OrH 

ooooooooooooo 

ooooooooooooo 


NrO'— i ro ^ i— ' (N m •/"} 'O 

>-> i>~> >^i >~. >~. >-. 

NNNNNNNNNNNN.. 


minmomcNcninmo^Hr-'CO 

o csj cn in o o cn o 

ooooooooooooo 


'sor-'sor-oo'or-ooo'* x x >> 

t^t^'OOOOOOO'.O'.ONO'. x >%!>% 

^ ^ . . 

P^fep^pHfetCHp^pHpHpHtlHpife 


OrOOt^O'^t^-OCNrO'— OO 
OOOCNOOOO’ — CS ’ — O (N 
OOOOOOOOOOOOO 





<t O r- ^TtCNOO'^trO^-t^-r^ 
in 0 ^- 00 <Nr<-) 0000 CNCN 
^‘oooooooooooo 


^j- i— 'CNrO^-m'sO'Or^'Ot^'OO'O 

■^t-inmirimiDvot^r^ooooooc^ 
0 \ 0 s 0 sG\G\ 0 s 0 \ 0 s 0 s 0 s 0 sG\ 0 ^ 
0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 0 \ 


(NOmoinin^j-cN^o^^HO 

P — ; — H 1 -H p P P p P O (N p 

ppooooop © © p © © 


O'vOsOsO'nO'nO'sOsO'nO'nOsOvO'nO'v 


00 (N 000 r- 0 ^H^H^H 0 (N 

p (N cn p (N co p (N O) p p p p 

000 ) 0000)000000 


i— • <n i— < <n m < <n m ^ ^ <n m ^t- 

oooooooooooooooooooooooooo 

oooooooooooooooooooooooooo 


ooi-nocNOcNinm^HOO^Hm 

PPP^H^H^HOPPPPP^H 

0 ) 0 ) 00000000 ) 0 ) 00 ) 


(Nro^J-iri'O^ot^-'— i ^ (N i— i<Nco 

i < N«Nromm 

t^r^r-t^t^r-r^oooooooooooo 

U*U*U*U*U*U*U*U*U-*U*U*U*U* 


^-HOCNmm^H^HOOOOcNO 

ooHTtiipt^yqp ^ p -h 
ooooooooooooo 

oo 11 III i I 


'O'sO'O'O'O'O'O'O'O'O'O'O'O 

V 0 ' 0 V > 0 V 0 V 0 ' 0 ' 0 V 0 ^ 0 '> 0 ' 0 V 0'0 


^tON^-rior-ooo^-ooo 

ooooooooooooo 

o> 1 in -xf 1 

CnIC *-! — ( Nrn ^ J - — — — <N — — <N 

— — <N rN — M rs 

PnfcpHpHfcpLiP^pHpHpHpLipHpH 



TOC Graphic 


